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SUMMARY 

Development of probabilistic structural analysis methodology for hot 
engine structures is a major activity at NASA Lewis Research Center. It con- 
sists of three program elements: (1) composite load spectra methods, (2) prob- 

abilistic structural analysis methods, and (3) probabilistic analysis and 
structural reliability application. Recent progress includes: (1) quantifica- 

tion of the effects of uncertainties for several variables on high pressure 
fuel turbopump (HPFT) turbine blade temperature, pressure, and torque of the 
space shuttle main engine (SSME), (2) the evaluation of the cumulative distri- 
bution function for various structural response variables based on assumed 
uncertainties in primitive structural variables, and (3) evaluation of the 
failure probability. Collectively, the results demonstrate that the structural 
durability of hot engine structural components can be effectively evaluated in 
a formal probabilistic/reliability framework. 


INTRODUCTION 

It is becoming increasingly evident that deterministic structural analysis 
methods will not be sufficient to properly design critical components in hot 
engine structures. These structural components are subjected to a variety of 
complex, and severe cyclic loading conditions including high temperatures and 
high temperature gradients. Most of these are quantifiable only as best engi- 
neering estimates. These complex loading conditions subject the material to 
coupled nonlinear behavior which depends on stress, temperature, and time. 
Coupled nonlinear material behavior is nonuniform, is very difficult to deter- 
mine experimentally, and perhaps impossible to describe deterministically. In 
addition, hot rotating structural components are relatively small. Fabrication 
tolerances on these components, which in essence are small thickness varia- 
tions, can have significant effects on the component structural response. Fab- 
rication tolerances by their very nature are statistical. Furthermore the 
attachment of components in the structural system generally differs by some 
indeterminant degree from that assumed for designing the component. In sum- 
mary, all four fundamental aspects: (1) loading conditions, (2) material 

behavior, (3) geometric configuration, and (4) supports - on which structural 
analyses are based, are of a random nature. The direct way to formally 
account for all these uncertain aspects is to develop probabilistic structural 
analysis methods where all participating variables are described by appropri- 
ate probability functions. 

NASA Lewis Research Center is currently developing probabilistic struc- 
tural analysis methods for select engine structural components (fig. 1). 
Briefly, the deterministic, three-dimensional, inelastic analysis methodology 
developed under the Hot Section Technology (HOST) and Research and Technology 
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Base Programs is being augmented to accommodate the complex probabilistic 
loading spectra, the thermovi scoplastic material behavior, and the material 
degradation associated with the environment of aerospace propulsion system 
structural components. The goal of the methodology is to address the problem 
depicted schematically in figure 2. 

The development of probabilistic structural analysis methodology consists 
of the following program elements (ref. 1): (1) composite load spectra models, 

(2) computational probabilistic structural analysis methods, and (3) probabil- 
istic analysis and structural reliability application. The development of the 
probabilistic structural analysis methodology is a joint program of NASA Lewis 
in-house and sponsored research. The objective of the proposed paper is to 
illustrate recent progress on the application of this methodology to determine 
the reliability of turbine blade components of rocket propulsion systems. Past 
progress of specific elements of the program are described in papers presented 
in several recent conferences (refs. 2 to 8). Recent activities have focused 
on extending the methods to include the combined uncertainties in several fac- 
tors on the structural response (fig. 3). The objective of the present paper 
is to briefly describe recent progress in three program elements: composite 

load spectra models, probabilistic finite element structural analysis, and 
probabilistic strength degradation modeling. Progress is described in terms 
of fundamental concepts, computer code development, and representative numeri- 
cal results. 


COMPOSITE LOAD SPECTRA 

The fundamental assumption for developing composite load spectra is that 
each individual load condition is the probabilistic time synthesis of four 
primitive parts: (1) steady state, (2) periodic, (3) random, and (4) spike. 

Each of these parts, except random, is described by a deterministic portion 
and a probabilistic perturbation about this deterministic portion as depicted 
schematically in figure 4. One justification for describing each loading 
condition in terms of primitive parts is that experts, over the years, have 
developed good judgments of the ranges of the perturbations about nominal 
(deterministic) conditions. The objective of the Composite Load Spectra pro- 
gram is to formalize the fundamental assumption in a computer code using: (1) 

available data from various rocket engines, (2) probability theory, and (3) a 
dedicated expert system. 

A schematic diagram of the composite load spectra (CLS) computer code is 
shown in figure 5. Representative results obtained for the perturbations of 
different engine factors on the SSME high pressure fuel turbopump turbine 
blade are shown in figures 6 and 7. Figure 6 depicts the nominal tempera- 
tures, while figure 7 indicates the temperature changes due to hot gas seal 
geometry and respective perturbations indicated in the figure caption. For 
example, the greatest temperature change due to gas seal geometry (fig. 7) is 
53.3 °F. The combined contributions of this and other factors (not shown 
here) is 87.6 °F which is in addition to the greatest nominal temperature of 
1860 °F. Although at first glance a change of 87.6 °F may seem insignificant, 
this is not the case, because at these levels, small temperature changes have 
dramatic effects on durability and cooling requirements. 
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Another representative example is shown in figure 8, where the effect of 
comparable factors on the torque of the high pressure oxidizer turbopump 
(HPOTP) are plotted as a function of time, for the start-up portion of a mis- 
sion cycle. These bounds are substantially wider at some times then at others 
Similar plots can be obtained for pressures or any other loading condition 
(ref. 5). The current CLS code permits the simultaneous perturbation of 47 
different factors for each different load condition. 


PROBABILISTIC FINITE ELEMENT STRUCTURAL ANALYSIS 

The fundamental assumption for developing probabilistic finite element 
structural analysis (PFEM) is that the uncertainties in each primitive struc- 
tural variable can be represented by an assumed probabilistic distribution. 
Primitive structural variables are those which are used to describe a struc- 
ture such as: (1) stiffness, (2) strength, (3) thickness and tolerance, (4) 

spatial location, (5) attachment, and (6) various nonlinear dependencies (tem- 
perature, stress, time, etc.). A schematic of the probabilistic distributions 
for some primitive variables is shown in figure 9. Subsequently, the uncer- 
tainties in the load conditions (described by the composite load spectra) and 
the uncertainties in the primitive structural variables are computationally 
simulated by performing probabilistic structural analyses to determine the ran 
dom structural response of a specific SSME structural component. The struc- 
tural response is generally described in terms of usual quantities such as 
displacement, frequencies, buckling loads, and structural fracture toughness. 
The integration is illustrated in figure 10. 

It is instructive to compare component development by the traditional 
engineering approach and component evaluation using PFEM. The parallelism is 
summarized in table I. The former approach relies on physical experiments for 
quantification and requires that the material, fabrication process, and test 
methods be available. The latter approach is entirely computational and 
requires the integration of available structural analysis methods with avail- 
able probability theory. The former approach has the advantage of demonstrat- 
ing a specific technology while the latter has the advantage of assessing 
undeveloped, but with high payoff potential, candidate technologies. In addi- 
tion, the former approach requires a large number of experiments to verify a 
design, while in the latter a design can be verified with strategically 
selected few experiments. 

PFEM has been formalized and integrated into a computer code identified 
as numerical evaluation of stochastic structures under stress (NESSUS). A 
schematic diagram of NESSUS is shown in figure 11. The user interacts with 
NESSUS through a dedicated expert system schematically shown in figure 12 
(ref. 8). Representative results obtained using NESSUS are shown in figure 13 
The distributions assumed for the primitive variables are listed in the table 
in figure 13. Both the individual and the combined effects of the primitive 
variable uncertainties on the combined stress (Von Mises) are shown in fig- 
ure 13 in terms of cumulative distribution functions (CDF). The information 
generated during the PFEM can be used to establish the confidence level. A 
sample result is shown in figure 14 for one blade location. 
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The curves in figure 14 can be used in a number of ways such as: (1) all 

the blades tested in the assumed conditions will have a mean combined stress 
between 57.9 and 62.1 ksi 90 percent of the time and (2) the combined stress 
in all the blades tested (under the assumed conditions) will range from 42 to 
83 ksi. This indicates that a wide scatter in the combined stress is proba- 
ble. Two implications that follow are: (1) assessing the durability/life of 

the blades using only material uncertainties will not be sufficient, and (2) 
obtaining wide scatter in measured stress/strain magnitudes does not neces- 
sarily indicate test procedure difficulties. 

All the NESSUS results presented herein were obtained using about 1000 
perturbations for each case studied. These perturbations are relatively small 
in number compared to direct Monte Carlo simulation which will normally require 
20 000 and greater simulations to generate the same curves (refs. 6 and 7). 

The reduced (but with comparable accuracy) number of simulations is a result 
of a NESSUS feature which utilizes the fast probability integration method 
(ref. 6) to select subsequent perturbations in a self-adaptive manner. 


PROBABILISTIC ANALYSIS FOR STRENGTH DEGRADATION 

The fundamental assumption for developing probabilistic analysis methods 
for strength degradation is that the uncertainties in the primitive variables 
affecting strength can be described by assumed distributions. Two different 
models are selected to demonstrate the concept. The models express the number 
of mechanical load cycles to failure. One of the models is based on linear 
elastic fracture mechanics and the other on a strength degradation model 
recently studied at NASA Lewis (ref. 9). The models with their respective 
primitive variables are summarized in table II. 

Both of these models are used to predict the number of cycles to failure 
in a material typical for SSME components. The input for the fracture mechan- 
ics model is summarized in table III. The CDF obtained from this input is 
shown in figure 15. The input for the strength degradation model is summa- 
rized in table IV and the corresponding CDF is shown in figure 16. Both CDF's 
exhibit wide ranges for the probable number of cycles to failure. The linear 
fracture mechanics model shows a mean of 10 000 cycles while the strength 
degradation model shows a mean of 10 million cycles. Based on this comparison 
the linear fracture mechanics model penalizes the material by three decades. 

It is important to note the differences between the two models: (1) the linear 

fracture mechanics model assumes the existence of a crack-like defect and then 
evaluates the number of cycles required to grow this defect to a critical size 
for imminent rapid propagation to fracture; (2) the strength degradation model 
does not presume the existence of defects and, therefore, includes both defect 
initiation and propagation resulting in greater number of cycles; (3) the lin- 
ear fracture mechanics model has five primitive variables while the strength 
degradation model has 13 (assuming that the greater the number of primitive 
variables in the model the more inclusive the representation of the physics in 
the model, then the strength degradation model will be more accurate); and 
(4) The linear fracture mechanics model requires determination of c, m, and 
a^ by specialized and often complex test methods while the strength degrada- 
tion model uses available room temperature material properties. The computer 
code to perform both of these simulations is described in reference 9. 
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CONCLUSIONS 


The development or probabilistic structural analysis methods for hot 
engine structures continues. Recent progress on its application to select 
SSME components includes: (1) the effects of the uncertainties of several fac- 

tors on blade temperatures, pressures, and torque, (2) the evaluation of the 
cumulative distribution function of structural response variables based on 
assumed uncertainties in the structural primitive variables, (3) evaluation of 
failure probability, and (4) life assessment in terms of cumulative distribu- 
tion function using linear fracture mechanics and strength degradation models. 
Three different computer codes are being developed in parallel: (1) Composite 

Load Spectra (CLS) for the probabilistic description of SSME component loads, 
(2) NESSUS, for the probabilistic structural analysis of select SSME structural 
components, and (3) a life durability code for the assessment of the fatigue 
cycles to failure of structural components in SSME mission environments. Col- 
lectively, the results obtained to date demonstrate that the structural dura- 
bility of SSME critical components can be evaluated using the probabilistic 
methodology under development. 
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TABLE I. - SOME HELPFUL PARALLELS 


Traditional component development 

Parameter identification 
Characteri zat i on - database 
Component fabri cation 
Component testing 
Component test results database 
Stati sti cal inferences 
Reliability/confidence level 

PSAM component evaluation 
Parameter i den t i f i ca t i on 
Probabilistic distributions 
Component mechanistic model 
Component analysis 
Component analysis results database 
Statistical inference 
Rel i abi 1 i ty/conf i dence level 

Bottom line 

Phys i cal 

Computati onal 


TABLE II. - PROBABILISTIC ANALYSIS FOR MATERIAL STRENGTH DEGRADATION AND FRACTURE M00ELS 
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TABLE III. - PRIMITIVE VARIABLE DISTRIBUTIONS FOR THE FATIGUE CRACK GROWTH MODEL 


Vari abl e 

Variable type 

Di stri bu- 
tion type 

Mean 

Standard deviation 

Material property, C 

Random 

Lognormal 

8. 66x10" 10 in. /cycle 

0.866x10"'° (10 percent) 

Material property, M 

Random 

Normal 

6.48xl0" 2 to./ 
cycl e/ksi i n . ' ^ 

3.24x1 0”^ (5 percent) 

Alternating range 




9 (10 percent) 

stress, A a 

Random 

Lognormal 

90 ksi 

Initial crack size, aj 

Random 

Lognormal 

118x1 0 -4 in. 

17.7x1 O'" 4 (15 percent) 

Final crack size, ap 

Deterministic 

N/A 

7.87xl0" 2 in. 

N/A 

Component/crack shape 

Determi ni st i c 

N/A 

1.0 

N/A 

parameter, Y 






TABLE IV. - PRIMITIVE VARIABLE DISTRIBUTIONS FOR THE STRENGTH REDUCTION MODEL 


Vari abl e 

Vari abl e 
Type 

Di stri but i on 
type 

Mean 

Standard deviation 

Melting temperature, Tp 

Random 

Normal 

2732 °F 

82.0 (3 percent) 

Ultimate tensile strength, Sp 



Lognormal 

130 ksi 

6.5 (5 percent) 

Log of final cycle, N M p 



Lognormal 

3 

.8 (10 percent ) 

Reference temperature, Tq 



Normal 

68 °F 

2.0 (3 percent) 

Residual compressive stress, oq 



Lognormal 

-2.9 ksi 

-0.145 (5 percent) 

Log of reference cycle, N^q 



Lognormal 

7 

0.7 (10 percent ) 

Reference fatigue strength, Sq 



Lognormal 

72.6 ksi 

3.6 (5 percent) 

Current temperature, T 



Normal 

1562 °F 

46.7 (3 percent) 

Current mean stress, o 



Lognormal 

21.8 ksi 

1.1 (5 percent ) 

Current fatigue strength, S 



Lognormal 

36.3 ksi 

1.8 (5 percent ) 

Temperature exponent, N 



Normal 

0.5 

0.15 (3 percent ) 

Stress exponent, M 



Normal 

0.5 

0.15 (3 percent) 

Cycle exponent, Q 



Normal 

0.5 

0.15 {3 percent ) 
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HIGH PRESSURE TURBOPUMP BLADE 



LOX POSTS 



MAIN COMBUSTION CHAMBER 


FIGURE 1. - SELECT SSME COMPONENTS 
INCLUDED IN THE PROBABLILISTIC 
STRUCTURAL ANALYSIS METHODS 
DEVELOPMENT PROGRAM. 
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ORIGINAL PAGE 

BLACK AND WHITE PHOTOGRAPH 


PROBLEM: 

[" NATURAL FREQUENCY FREQUENCY ] 

prob |^ 0F BLAD £ / o ^ limit, n L J = ? 



CONFIDENCE 

• DIFFICULT ANALYSIS 
PROBLEM 

• APPROXIMATE METHOD 
DEFINED 

• USES NESSUS DATA BASE 

• GIVES INTERVAL ESTIMATES 

- PROBABILITY 

- VARIABILITY 


FIGURE 2. - PROBABILISTIC STRUCTURAL ANALYSIS DEFINITION. 


PLANNED ACCOMPLISHMENT 


STATUS 


DEVELOP PROBABILISTIC STRUCTURAL 
ANALYSIS METHODS/COMPUTER CODES 
FOR QUANTIFYING THE INDIVIDUAL 
AND COMBINED EFFECTS OF UNCER- 
TAINTIES ON THE STRUCTURAL DUR- 
ABILITY OF SSME COMPONENTS 


COMBINED PROBABILITY OF OCCURRENCE 



BLADE STRESS/ KSI 


• FORMULATED PROBABILISTIC FINITE 
ELEMENT METHODS 

• DEVELOPED SECOND GENERATION 
COMPUTER CODE - OPERATIONAL AT 
LEWIS. UNCERTAINTIES INCLUDED: 

- LOADS 

- GEOMETRY 

- MATERIAL 

- SPATIAL POSITION 

- ATTACHMENT 

SIGNIFICANCE 

THE COMBINED PROBABILITY OF MAG- 
NITUDE OF STRUCTURAL DURABILITY 
VARIABLES CAN BE DETERMINED 


FIGURE 3. - PROBABILISTIC STRUCTURAL ANALYSIS METHODS (PSAM). 
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TIME CD-0S-321&4 

FIGURE H. - EACH LOAD CONDITION IS THE PROBABILISTIC TIME 
SYNTHESIS OF FOUR PRIMITIVE PARTS: COMPOSITE LOAD SPECTRA. 


COMPOSITE LOAD SPECTRA EXPERT SYSTEM 


LOAD EXPERT SYSTEM DESIGN PHILOSOPHY 


EXPERT 

SYSTEM 



• RULE-BASE PRODUCTION SYSTEM 

• IF-THEN RULES 

• SIMPLE INFERENCE SCHEME 

• INFERENCE NET (DECISION TREE) 

• SOPHISTICATED PROBABILISTIC METHODS 

• DISCRETE PROBABILITY DISTRIBUTION 

• MONTE CARLO 

• BARRIER CROSSING 


• POWERFUL KNOWLEDGE BASE 

• INFLUENCE COEFFICIENTS 

• SCALING COEFFICIENTS 

• NOMINAL ENGINE COEFFICIENTS 

• DUTY CYCLE LOAD PROFILES 

• ENGINE CONFIGURATION AND GEOMETRY DATA 

• RAW ENGINE FLIGHT AND TEST DATA 

FIGURE 5. - COMPOSITE LOAD SPECTRA SIMULATION USING EXPERT SYSTEMS. 
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ORIGINAL PAGE 

ACK AND A/HITE PHOT€KiRA£H 


STEADY STATE TEMPERATURE 


FIGURE 6. - COMPOSITE LOAD SPECTRA PREDICTED STEADY STATE TEMPERA 


TEMPERATURE PERTURBATION DATA DUE TO HOT GAS SEAL GEOMETRY 
DELTA TEMPERATURES 
ONE STANDARD DEUIATION OF .060 


FIGURE 7. - TEMPERATURE PERTURBATION DATA DUE TO HOT GAS SEAL 
GEOMETRY DELTA TEMPERATURES. ONE STANDARD DEVIATION. 



TORQUE, ft-lb, THOUSANDS 


1*2 I — 







ATTACHMENT 

MAGNITUDE 



DEPENDENCE 


CO-88-32183 


FIGURE 9. - PROBABILISTIC STRUCTURAL ANALYSIS: THE UNCER- 
TAINTIES IN EACH STRUCTURAL PRIMITIVE VARIABLE CAN BE 
DESCRIBED BY ASSUMED CORRESPONDING PROBABILISTIC DIS- 
TRIBUTIONS. 


r NODAL LINE 



RESPONSE VARIABLES (u, w, a, t, G) ARE 
DETERMINED FROM THE RESULTS OF THE 

MULTIPLE FINITE ELEMENT ANALYSES CD . e8 . 32187 

FIGURE 10. - UNCERTAINTIES IN THE LOAD CONDITIONS AND STRUCTURAL 
PRIMITIVE VARIABLES ARE INTEGRATED BY PERFORMING MULTIPLE FINITE 
ELEMENT ANALYSES. 


12 


FINITE ELEMENT LIBRARY 


PERTURBATION VARIABLES: 


DATA FLOW THROUGH THE 
MAIN PROCESSORS IN NESSUS 


a & & 




PLANE STRESS PLANE STRAIN AXISYMMETRIC 



3D BRICK 3D SHELL LINEAR BEAM 


• ALL ELEMENTS INTEGRATED NUMERICALLY 

• LINEAR LAGRANGIAN INTERPOLATING 

FUNCTIONS 

• SELECTIVELY REDUCED INTEGRATION 


(I) TOPOLOGY 

NODAL COORDINATES 
NODAL SHELL THICKNESS 
NODAL SHELL OR BEAM NORMALS 

(II) GEOMETRIC PROPERTY 

THICKNESS OF PLAIN STRESS 
ELEMENTS 

ELASTIC BEAM SECTION PROPERTIES 

(III) MATERIAL PROPERTY 

ELASTICITY MODULUS 
POISSON'S RATIO 
THERMAL EXPANSION COEFFICIENT 
MATERIAL DENSITY 

(IV) BOUNDARY CONDITIONS 

BASE SPRING STIFFNESSES 

(V) LOADS 



NODAL FORCE VECTORS 
ELEMENT PRESSURES AND EDGE 
TRACTIONS 

NODAL TEMPERATURES 


FIGURE 11. - PROBABILISTIC FINITE ELEMENT COMPUTER CODE NESSUS. 
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RANDOM VARIABLE DISTRIBUTION TABLE 



i x 

S-L SSME 
BLADE 


CASE 

RANDOM 

VARIABLE, 

RV 

a COV OF RV 

1 

X - COOR (N) b 

0.005 (STD) C 

2 

Y - COOR (N) 

.005 (STD) 

3 

Z - COOR (N) 

.005 (STD) 

4 

T - TEMP (N) 

60 °F (STD) 

5 

| 

E - YOUNG'S 
MODULUS (N) 

0.10 

6 1 

i 

L - LATERAL 
PRESSURE (N) 

0.20 

7 

F - CONC. 
LOAD CN) 

0.20 

8 

9 j 

R - ROTATION 
SPEED (N) 

C - COMBINE 

0.02 


CDF OF VON MISES STRESS DUE TO SINGLE/ 
COMBINED RANDOM VARIABLE PERTURBATION 


Z COORDINATE 



1 .0 j — 

.5 - 
0 — 


LATERAL 

PRESSURE 


YOUNG'S 

MODULUS 


ROTATION 

SPEED 


CONCENTRATED 

LOAD 


1.0 

.5 


X COORDINATE 


TEMPERATURE 


r r r 

i 


COMBINED CASE 


0 60 120 0 60 120 0 60 120 
COMBINED STRESS AT ROOT, KSI 


a COV: COEFFICIENT OF VARIATION ( STD/MEAN) 

b (N) NORMAL DISTRIBUTION 
C STD: STANDARD DEVIATION 


FIGURE 13. - PROBABILISTIC STRUCTURAL ANALYSIS USING NESSUS. 


o 

s 



COMBINED STRESS 


FIGURE 14. - 90% CONFIDENCE BOUND RESULT FROM 
UNCERTAINTY IN THE INPUT PARAMETERS. 
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PROBABILITY FOR FRACTURE STRESS 



FIGURE 15. - SIMULATED FAILURE PROBABILITY USING A FATIGUE 
CRACK GROWTH MODEL. 



FIGURE 16. - SIMULATED FAILURE PROBABILITY USING 
A STRENGTH REDUCTION MODEL. 
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